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Wavelet  Based  Methods 


Wavelet  image  representations 


The  ability  of  wavelets  to  localize  position  and  frequency  is  very  useful  in  representing 
images. 

This  allows  us  to  isolate  a  region  of  the  image  or  to  examine  the  image  on  different  scales. 

Decomposing  an  image  in  a  wavelet  basis  tends  to  concentrate  most  of  the  energy  of  the 
image  into  a  few  (large)  coefficients.  This  is  referred  to  as  energy  compaction. 

One  reason  for  this  is  that  images  tend  to  have  large  smooth  regions.  In  these  regions  the 
wavelet  coefficients  decay  rapidly. 

This  also  means  that  the  places  where  the  wavelet  coefficients  fail  to  decay  rapidly  (usually) 
correspond  to  edges  or  texture. 
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JPEG2000 


The  new  JPEG2000  standard  uses  wavelet  coding  for  image  compression  (both  lossless  and 
lossy). 

In  the  lossless  mode,  an  integer-to-integer  wavelet  transform  is  used.  This  mode  can  also  be 
lossy  if  we  cut  the  progressive  bitstream  at  some  point. 

The  coefficients  are  quantized  and  then  each  bitplane  is  coded  independently. 

When  a  particular  bit  (in  a  bitplane)  is  encoded  depends  on  its  “significance”.  Thus,  even 
within  a  bitplane  the  stream  can  be  progressive. 

A  context-dependent  binary  arithmetic  coder  is  used  for  entropy  coding. 

The  progressive  local  improvement  data  is  collected  in  “packets”  which  form  “layers”  (a 
global  quality  improvement). 

The  scale  and  resolution  capability  of  wavelets  is  essential  for  this  progressive  encoding. 
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Separable  decomposition 


.  -  p.6/7; 


Non-separable  wavelets 


The  scaling  function  for  a  non-separable  wavelet  satisfies  a  matrix  dilation  equation  of  the 
form 

4>(x)  =  Ci<j>(Ax  —  i ) 

where  A  is  a  dilation  matrix  (expansive  integer  matrix). 

The  simplest  examples  of  non-separable  wavelets  are  Haar  wavelets  based  on  self-affine 
tilings  of  JR2. 
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Other  self-affine  tilings 


Non-separable  subsampling 
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Smooth  non-separable  wavelets 


The  construction  of  continuous  or  smooth  non-separable  wavelets  is  considerably  more 
difficult  than  Haar  wavelets. 

There  are  several  obstructions  to  using  the  one-dimensional  techniques  to  construct  such 
wavelets. 

At  the  present  time,  not  many  families  of  smooth  non-separable  wavelets  are  known.  It  is 
known  that  arbitrarily  smooth  non-separable  compactly  supported  wavelets  exist. 


Truncation:  separable  vs.  non-separable 


Non-separable  truncations  6,8,10,12  levels 
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Truncation:  separable  vs.  non-separable 


Separable  truncations  2, 3, 4,5  levels 


Non-separable  truncations  4,6,8,10  levels 
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Wavelet  Denoising 


Suppose  we  have  an  image  X  that  is  corrupted  with  noise  N  so  we  observe  Y  =  X  +  N. 
We  want  to  find  an  estimate  X  for  X. 

The  energy  compaction  property  of  wavelets  means  that  Y  will  have  a  few  large  coefficients 
and  many  small  coefficients. 

The  basic  idea  in  wavelet  denoising  is  to  find  these  small  coefficients  and  to  shrink  them 
(either  zero  them  or  shrink  them  towards  zero),  thus  reducing  the  noise. 
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Hard  vs.  Soft  Thresholding 


There  are  several  different  methods  for  setting  the  threshold.  One  basic  distinction  is  hard 
versus  soft  thresholding. 


Soft  thresholding  has  the  advantage  of  being  continuous  (so  as  not  to  introduce  artificial 
discontinuities). 

Soft  thresholding  has  also  been  shown  to  achieve  a  near-optimal  minimax  rate  over  several 
Besov  spaces. 

Soft  thresholding  also  has  some  other  theoretical  advantages. 


Denoising  (cont) 


BayesShrink:  the  threshold  is  <r2 /crX-  This  threshold  is  derived  by  assuming  that  the 
distribution  of  the  wavelet  coefficients  in  any  scale  is  given  by  a  Generalized  Gaussian 
Distribution  (GGD).  The  problem  is  to  estimate  the  GGD  parameters  for  each  subband. 
These  parameters  are  estimated  for  each  subband  from  the  noisy  image  data. 

SureShrink:  the  threshold  is  computed  as  either  a  SURE  threshold  (derived  from  minimizing 
Stein’s  unbiased  risk  estimate)  or  the  Universal  threshold  (a  a/2  log(M),  which  is 
asymptotically  optimal  in  the  minimax  sense). 

The  SURE  threshold  can  be  viewed  as  another  (and  data-driven)  way  to  minimize  the  Bayes 
risk  given  a  GGD  prior  on  the  coefficients. 

NormalShrink:  the  threshold  is  (3cr2 /cry ,  where  (3  =  ^/\og(Lk/ J)  ( Lk  is  the  length  of  the  k\h 
subband  and  J  is  the  number  of  subbands),  a  is  the  estimated  noise  variance  (as  before) 
and  cry  is  the  variance  of  the  particular  subband  of  the  noisy  signal. 


.-p.15/7: 


Cycle  Spinning 


The  discrete  wavelet  transform  has  a  lattice  of  preferred  locations.  These  locations  might  not 
be  the  correct  ones  for  denoising  a  particular  image. 

One  way  to  deal  with  this  problem  is  via  cycle  spinning-  denoising  shifted  versions  of  the 
image  and  then  averaging  the  various  shifts. 


The  idea  is  that  the  noise  at  different  locations  is  independent,  so  the  averaging  will  tend  to 
smooth  out  the  noise. 
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IFS  Based  Methods 


IFS  Fractal  Basics 


We  wish  to  formalize  the  notion  of  a  set  being  self-similar. 
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Consider  the  three  affine  maps 

/  \  ( •T'  V  \  t  ^  ,X  +  1  y  X  y  +  1 

wo{x,y)  =  W!  (x,y)  =  w2(x,y)  =  ) 

We  see  that  if  5  is  the  Sierpinski  gasket  (above)  then 

S  =  wq(S)  U  wi(S)  U  w2(S). 


In  fact,  S  is  the  unique  such  set. 
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Definition  of  IFS 


Let  (X,  d)  be  a  complete  metric  space. 

An  IFS  on  X  is  a  finite  set  of  contractive  self-maps  {wi  :  X  ^  X  :  i  =  0,. .  .N}. 
The  attractor  of  the  IFS  is  a  non-empty  compact  set  A  c  X  so  that  A  =  (J*  Wi(A). 
We  endow  the  space 

H(X)  =  {K  c  X  :  Q)  ^  K  \s  compact} 
with  the  Hausdorff  metric  dh  making  (' U(X),dh )  a  complete  metric  space. 

For  our  purposes,  we  usually  take  X  =  [0, 1] 2 . 


.-p.19/7: 


The  map  W  :  7i(X)  'U(X)  defined  by 


W(B)=\Jwi(B) 

i 


is  a  contraction  if  {wi}  is  an  IFS. 

Thus 

Bn  :=  Wn(B)  A 


for  any  non-empty  compact  set  B. 
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Collage  Theorem 


Theorem:  (Collage  Theorem)  Let  T  :  X  X  be  a  contractive  map  with  contractivity 
factor  s  <  1  and  xT  be  its  fixed  point.  Then 


d(y,xT)  < 


d(y,T(y)) 
1  —  s 


In  the  context  of  an  IFS  {wi},  this  says  that  if  a  set  B  almost  self-tiles  under  {wi}  then  B 
“close”  to  the  attractor  of  the  IFS. 

This  is  used  in  the  inverse  problem:  Given  a  set  A,  find  an  IFS  whose  fixed  point  is  A. 
Often  the  quantity  d(y,T(y))  is  referred  to  as  the  collage  distance. 


Collage  Theorem  (cont) 


A  collage  of  a  leaf  with  copies  of  itself. 


The  attractor  of  the  resulting  IFS. 
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IFS  on  Function  Spaces  (IFSM) 


We  consider  an  image  as  a  function  /  :  X  1R,  so  we  need  an  IFS  operator  on  functions. 
Let  4>i  :  1R  1R  be  maps.  We  usually  want  them  to  be  Lipschitz. 

Define  T(f)(x)  =  ^  (/(^_1(x)))  so  T  maPs  functions  on  X  to  functions  on  X. 

i 

Under  suitable  conditions  on  the  wi’s  (contractivity)  and  on  the  </>*’ s  we  know  that  T  will  map 

LP(X)  to  lp(X). 

Notice  that  the  w^s  perform  a  mapping  on  the  “base  space”  X  and  the  </>*’s  map  on  the 
“vertical”  space  R  The  ^’s  are  usually  called  the  grey-scale  maps. 

Most  often  we  take 

Wi(x)  =  dix  +  Si  and  <j>i(t)  =  <xit  + fa. 
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IFSM  example 


There  is  one  parent  and  four  children. 


a 

0.5 

0.5 

0.4 

0.5 

p 

80 

40 

132 

0 

The  parameters  are 


Another  example 


This  is  the  same  type  of  IFSM  operation,  except  with  a  twist. 


Here  we  have  allowed  non-linear  maps  wi. 


Fractal  Image  Compression 


Given  a  picture  (a  function  /),  the  basic  idea  is  to  try  to  find  an  IFSM  operator  T  whose  fixed 
point  xT  closely  matches  /. 

If  we  can  find  a  T  so  that  d(T(f),f)  is  small,  then  we  know  (by  the  Collage  Theorem)  that 
d(f,xT)  will  be  small  as  well. 

The  standard  fractal  compression  method  is  the  Block-Encoding  method. 


The  maps  u>i  are  found  by  the  searching  while  the  a*  and  ft  are  (usually)  calculated  to 
minimize  mean  square  error. 


Fractal  Compression  as  Self-VQ 


Vector  Quantization  (VQ)  is  one  form  of  data  compression  for  multi-dimensional  data. 

In  basic  VQ  you  map  each  data  vector  into  an  appropriate  (usually  closest)  code  vector  from 
a  codebook. 

This  allows  you  to  quantize  vector  data  which  might  have  dependencies  between  the 
components. 


IFS  block-coding  is  like  a  Self-VQ  in  that  the  image  acts  as  its  own  codebook. 


Why  are  we  interested  in  compression? 


Compression  is  approximation  -  we  attempt  to  approximate  the  image  by  the  fixed  point  of 
an  IFSM. 

In  doing  this,  we  analyze  the  “fractal”  nature  of  the  image. 

The  starting  point  of  the  IFS-based  methods  for  image  processing  is  this  fractal  analysis. 


.  -  p.28/7: 


Local  IFS 


This  type  of  IFS  can  be  formalized  as  a  Local  IFS. 


-  p.29/7: 


Basic  Fractal  Compression:  example 


Reconstruction:  one,  two,  three,  four,  ten  iterations  and  the  original  image. 


.  -  p.30/71 


Extensions  to  the  basic  algorithm 


Many  modifications  and  extensions  to  the  basic  algorithm  have  been  proposed. 

These  include  adaptive  partitions  (both  top  down  and  bottom  up),  more  complicated 
grey-level  maps  (the  fa's)  and  merging  fractal  and  wavelet  or  other  transform  techniques. 

A  significant  problem  is  the  problem  of  searching  for  the  best  match.  This  has  also  been  a 
large  focus  of  the  research  in  this  area. 

Methods  to  speed  up  the  search  include  classifying  the  blocks,  limiting  to  a  local  search, 
using  transforms  or  convolution,  and  doing  the  search  in  a  hierarchical  fashion  (multi-scale). 
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Postprocessing 


One  problem  with  the  iteration  scheme  is  that  edges  (and  noise)  can  get  replicated  to  finer 
scales. 


One  can  either  apply  a  smoothing  filter  to  the  end  result  or  (better)  apply  some  smoothing 
operation  for  each  iteration. 


.  -  p.32/7: 


Varying  the  partition 


Many  different  types  of  partitions  have  been  tried,  including  triangular  partitions,  rectangular 
partitions  and  more  general  partitions. 

Most  of  these  partitioning  schemes  are  also  adaptive. 

The  benefit  is  the  flexibility  to  match  image  features. 

The  cost  is  both  the  extra  algorithmic  complexity  and  the  extra  storage  necessary  for 
encoding  the  partition. 
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Examples  of  three  partitions 


Quadtree  partition 


This  is  the  simplest  adaptive  partitioning  scheme. 

Start  with  large(ish)  child  blocks  and  try  to  find  a  match,  subdividing  into  four  smaller  blocks  if 
some  error  tolerance  is  not  met. 

The  tree  structure  allows  simpler  storage  of  the  partition. 
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HV  partitions  (Horizontal,  Vertical) 


There  are  several  split  criteria. 

We  could  try  to  find  an  edge  and  split  on  the  edge. 

We  could  split  to  minimize  the  error  (often  mean-square). 


.  -  p.36/7: 


Triangular  partition 


Triangular  partitions  are  very  flexible  but  complex. 
The  partition  requires  a  lot  of  information  to  encode. 
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Region  growing 


In  this  approach,  we  start  with  a  partition  into  small  atomic  regions,  usually  blocks.  These 
blocks  function  as  the  child  regions. 

We  selectively  merge  two  regions  according  to  an  error  criteria. 

We  must  also  merge  the  corresponding  parent  regions. 

Thus,  we  keep  a  list  of  active  child  and  parent  regions  and  possible  candidates  for  merging. 
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Region  growing  example 


Region  growing  partition  and  quadtree  partition 


Region  growing  example 


Variance  based  splitting  criteria  Two  collage  based  splitting  criteria 
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Beyond  collage  coding 


The  IFS  image  approximation  usually  uses  the  collage  distance  as  a  measure  of  error. 


d(y,xT)  < 


d(y,T(y)) 
1  —  s 


This  provides  a  sub-optimal  solution,  since  we  really  want  to  minimize  d(y,  xT)- 
This  is  done  to  simplify  the  computational  problem. 

There  are  various  ways  of  trying  to  improve  on  this. 

These  usually  involve  starting  with  the  collage  solution  and  doing  some  type  of  “local”  search 


Gradient  descent  optimizing 


Unfortunately,  gradient  descent  doesn’t  really  pay  off  in  terms  of  image  quality  and  it  is  very 
computationally  expensive 
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F ractal  zoom 


.  -  p.43/7: 


F ractal  Denoising 


Most  approaches  to  denoising  involve  some  type  of  local  smoothing,  with  the  assumption 
that  the  noise  is  additive  and  zero  mean. 

Using  fractal  compression  as  a  basis  for  denoising  is  natural  since  fractal  compression 
involves  some  local  smoothing. 

This  smoothing  arises  from  two  sources:  the  mean-square  fit  in  the  region  matching  and  the 
contractivity  in  the  grey-scale  maps,  </>*. 

We  can  also  try  to  estimate  the  noise-free  fractal  parameters  and  use  these  in  the  fractal 
denoising. 

Any  of  the  fractal  compression  schemes  can  be  used  for  denoising. 

For  this  application,  we  don’t  really  care  about  compression  so  approximation  schemes 
which  expand  are  ok. 
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Example:  Fractal  Denoising 


Original 


Standard  scheme,  cycle  spinning  with  16  shifts 


Noisy 


Quadtree  scheme,  cycle  spinning  with  16  shifts 


Mixing  IFS  and  wavelets  (IFSW) 


The  key  observation  in  mixing  IFS  and  wavelets  is  the  fact  that  wavelets  have  a  natural 
notion  of  scale. 

We  expand  our  signal  in  some  wavelet  basis  and  then  try  to  do  an  IFS  type  compression  on 
the  wavelet  coefficients. 


For  simplicity,  we  restrict  our  discussion  to  one  dimension.  The  extension  to  two  dimensions 
is  simple. 

The  type  of  transform  considered  is  illustrated  as 


Basically  we  are  trying  to  extrapolate  the  fine  resolution  wavelet  coefficients  from  the  coarse 
resolution  wavelet  coefficients. 

Usually  one  will  store  all  the  expansion  coefficients  up  to  some  level  and  then  try  to  use  an 
IFS  to  predict  the  higher  frequency  ones. 


l 
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One  of  the  main  motivations  for  this  is  that  the  wavelets  overlap,  thus  the  “blocks”  in  the 
“wavelet  coefficient  space”  are  overlapping  in  “picture  space”. 

However,  since  the  wavelets  are  orthogonal,  the  different  “blocks”  still  represent  independent 
information. 

This  helps  to  eliminate  some  of  the  blockiness  from  the  fractal  block  algorithm. 

Wavelets  are  localized  in  both  frequency  and  space,  thus  different  “blocks”  correspond  to 
different  parts  (spatially)  of  the  picture. 

If  we  use  the  separable  Haar  basis,  IFSW  compression  is  exactly  the  same  as  IFS 
compression. 

For  other  wavelets,  the  relationship  is  more  complicated. 

One  could  also  try  to  do  an  IFS-type  compression  scheme  in  other  transform  spaces  (or 
using  other  orthonormal  basis). 
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Comparison :  F ractal/F ractal- Wavelet 


Standard  IFSM  scheme 


IFSW  scheme 
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Extensions  to  IFSW 


There  are  several  natural  ways  to  generalize  the  basic  IFSW  transform. 
The  2D  wavelet  coefficients  are 
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We  can  picture  a  general  2D  IFSW  as 
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Introduction  to  CBIR 


Content  Based  Image  Retrieval  is  an  application  where  you  query  an  image  database  and 
the  database  returns  the  N  most  “visually  similar”  images. 

Of  course,  the  major  difficulty  is  defining  what  “similar”  means  -  this  is  an  ill-posed  problem. 

The  other  problem  is  one  of  speed  -  the  query  should  take  only  a  couple  of  minutes  even 
with  databases  of  100  000  images  or  more.  This  speed  restriction  severely  constrains  the 
possible  image  analysis. 

The  usual  method  is  to  calculate  “image  descriptors”  -  these  are  typically  vectors  of  floating 
point  numbers. 

These  same  descriptors  are  computed  for  the  query  image  and  compared  to  those  for  the 
images  in  the  database.  A  score  is  calculated  for  each  image  and  the  top  N  images  are 
returned. 

The  mathematical  challenge  is  to  discover  descriptors  and  metrics  to  compare  the 
descriptors. 
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Example:  Color  Histogram 


As  an  example,  a  simple  image  descriptor  is  a  color  histogram. 

To  compute  this  descriptor,  choose  some  number,  say  24,  of  representative  colors. 

Map  all  the  colors  in  an  image  to  one  of  these  representative  colors. 

Scan  through  all  the  pixels  in  the  image  and  record  how  many  times  each  representative 
color  is  used. 

We  now  must  compare  two  color  histograms  to  compute  a  distance  between  two  images. 
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We  could  use  an  lp  metric,  but  this  has  disadvantages.  For  some  purposes,  an 
Earth-Movers  distance  is  better. 


All  lp  metrics  would  measure  these  distances  as  the  same,  but  in  color  space  blue  might  be 
closer  to  purple  than  yellow. 

Thus,  we  want  a  metric  that  measures  both  the  “horizontal”  difference  as  well  “vertical” 
difference. 
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F  ractal  Descriptors 


During  fractal  compression,  a  lot  of  data  about  the  image  is  generated. 

This  information  can  be  used  to  try  to  classify  the  image. 

A  very  simple  example  of  this  would  be  an  overall  collage  error.  That  is,  once  the  final  fractal 
“code”  T  has  been  found  for  the  image  y,  compute  d(y,  T(y)). 

This  gives  a  measure  of  how  self-similar  the  image  is  (at  least  under  this  class  of  IFS). 

One  could  also  compute  a  histogram  of  “local  collage  errors”-  that  is,  the  matchings  for  each 
parent-child  block  matching. 

Other  possibilities  include  statistics  of  relative  location  of  parent-child  matchings,  statistics  of 
a  and  (3  parameters,  etc. 
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F ractal  Descriptors  (cont) 


The  meaning  of  these  fractal  descriptors  is  not  completely  clear,  but  they  do  provide  one  way 
of  measuring  similarity. 

Which  metric  to  use  for  which  descriptor  is  also  not  at  all  clear. 

Edge  information  is  partially  stored  in  the  parent-child  matchings. 

The  (3’s  partially  encode  the  “global  look”  of  the  image.  They  represent  something  like  the 
average  value  of  the  image  over  the  child  block. 
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Meaning  of  IFS  parameters 


Fixed  a,  (3,  index, 

random  a,  (3 ,  index,  reconstructed. 
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IFS  parameters  (cont) 


Fixed  (a,  (3),  (a, index),  (/3, index), 

random  (a,/3),  (a, index),  (/3, index),  reconstructed 


Example:  Image  Matching 
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Example:  Image  Matching 
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Example:  Image  Matching 
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Example:  Image  Matching 
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Finding  subimages 
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Multifractal  Based  Methods 


Basics 


“The  multifractal  analysis  of  images  consists  in  defining  a  sequence  of  capacities  on  the 
image,  computing  its  multifractal  spectrum,  and  classifying  each  point  according  to  the 
corresponding  value  of  (a,  /(a)),  both  in  a  geometric  and  a  probabilistic  fashion.”  -  Jacques 
Levy-Vehel 


A  point  on  an  edge  (piecewise  smooth  curve)  in  the  image  belongs  to  a  set  Ea  whose 
dimension  is  1. 

A  point  in  a  homogeneous  region  has  f(a)  =  2. 

The  idea  is  that  edge  points  should  be  relatively  rare  points  -  too  many  edge  points  in  a 
region  should  be  classified  as  a  texture,  not  a  collection  of  edges. 

This  allows  us  to  analyze  the  image  in  several  ways. 
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Some  Examples  of  Capacities 


Sum  measure:  For  a  set  Q  we  simply  sum  over  the  image  grey  level  values  in  the  set. 

Max  capacities:  We  can  define  these  at  different  resolutions.  Basically,  we  compute  the 
maximum  of  the  pixels  in  a  region  Q  or  the  maximum  block  (at  some  resolution)  which 
intersects  Q. 

Iso  capacities:  We  can  again  define  these  at  different  resolutions.  These  are  defined  to  be 
the  size  of  the  largest  level  set  in  Q  at  some  resolution. 

The  max  capacities  depend  only  on  the  grey  level  values  while  the  iso  capacities  depend 
only  on  the  grey  level  distribution.  Thus,  these  give  different  information. 
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Edge  Detection:  Example 


Original  Image,  Smooth  Edges  obtained  with  max  capacity,  Irregular  Edges  obtained  with 
max  capacity. 

We  define  a  point  to  be  on  a  “smooth”  edge  if  it  belongs  to  a  set  Ea  whose  dimension  is  1. 

We  define  a  point  to  be  on  an  “irregular”  edge  if  it  belongs  to  some  set  Ea  whose  dimension 
is  between  two  values,  say  1.1  and  1.5. 

Edge  points  are  defined  using  both  local  information  (a(x))  and  global  information  (/(«)). 
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Multifractal  Denoising 


The  basic  idea  in  using  multifractal  analysis  in  image  denoising  is  that  we  want  to  increase 
the  Holder  function  a  of  the  image. 

We  want  to  do  this  in  such  a  way  as  to  preserve  the  important  features  of  the  image, 
particularly  edges. 


The  denoised  signal  X  is  thus  viewed  as  the  image  closest  to  Y  which  has  the  desired 
Holder  function, 

Usually  we  take  ax  =  oly  +  S,  where  S  >  0  is  some  function  (so  that  X  everywhere  has 
greater  smoothness  than  Y). 
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An  immediate  problem  arrises:  if  the  noise  is  additive  white  noise  then  ay  =  -1/2  a.s. 
everywhere.  Thus  there  is  no  way  to  extract  information  about  ax  from  aY- 

The  solution  to  this  is  that  we  have  only  a  limited  number  of  resolutions.  We  craft  our 
measures  of  regularity  to  agree  with  our  intuition  at  these  resolutions. 

Given  this,  we  still  have  to  be  able  analyze  and  prescribe  the  local  Holder  regularity. 

One  solution  to  both  problems  involves  wavelets.  There  are  simple  estimates  on  the  wavelet 
coefficients  which  guarantee  a  certain  local  Holder  regularity: 


\did\  <  C2~^a+1/2\ 


Thus,  given  the  wavelet  coefficents  dij  of  Y  we  modify  them  to  obtain  the  coefficients  cij  of 
X,  hopefully  with  the  desired  local  smoothness  and  still  close  to  Y. 
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Prescribing  the  Holder  Function 


Each  wavelet  coefficient  at  a  “large”  resolution  influences  many  points.  For  this  reason,  the 
minimization  problem  is  a  global  one  -  we  must  solve  simultaneously  for  all  the  wavelet 
coefficients  as  one  large  optimization  problem. 

We  obtain  a  system  of  constraint  equations  which  is  linear  in  the  logarithm  of  the  desired 
wavelet  coefficients.  Our  objective  function  is  the  risk  E(\\X  -  Y||). 

To  simplify  the  situation,  we  could  assume  that  Cij  =  Bjdij  with  0  <  Bj  <  1.  Thus,  we 
shrink  the  coefficients  of  Y  to  get  those  of  X.  Notice  we  allow  ourselves  to  shrink  each 
resolution  by  a  different  amount. 
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Original  SAR  image,  Kuan  median  filtering, 


Original  SAR  image  and  local  regularity  based  denoised  version. 
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Resources 


-•  Fractal  Image  Compression,  Yuval  Fisher,  editor,  Springer- Verlag,  1995. 
-•  The  Konstantz  Paper  Collection, 

www. inf.uni-konstanz .de/cgip/ f ractal2 /html / index . html 

-•  Waterloo  Fractal  Compression  Project,  links . uwaterloo .  ca 
&  Jacques  Levy-Vehel’s  web  page  at  f ractaies  .  inria .  f r. 

-•  “A  wavelet  tour  of  signal  processing”,  Stephane  Mallat,  Academic  Press. 
M  “Introduction  to  the  multifractal  analysis  of  images”,  Jacques  Levy-Vehel. 
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